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Abstract 

We present a polymer quantization of spherically symmetric Einstein gravity in 
which the polymerized variable is the area of the Einstein- Rosen wormhole throat. 
In the classical polymer theory, the singularity is replaced by a bounce at a radius 
that depends on the polymerization scale. In the polymer quantum theory, we 
show numerically that the area spectrum is evenly-spaced and in agreement with 
a Bohr-Sommerfeld semiclassical estimate, and this spectrum is not qualitatively 
sensitive to issues of factor ordering or boundary conditions except in the lowest 
few eigenvalues. In the limit of small polymerization scale we recover, within the 
numerical accuracy, the area spectrum obtained from a Schrodinger quantization 
of the wormhole throat dynamics. The prospects of recovering from the polymer 
throat theory a full quantum-corrected spacetime are discussed. 
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1 Introduction 



Polymer representation of quantum mechanics [H Ej has recently attracted con- 
siderable interest as a way to model loop quantum gravity in a technically simple 
setting. It has been applied to quantum mechanical systems including the har- 
monic oscillator [Tj, Coulomb potential [3] and 1/r 2 potential [I]. With these 
singular potentials it avoids the classical singularity in a way that is distinct from 
that in which the singularity is absent in the usual Schrodinger quantization. 
When applied to gravitational mini-superspace models, polymer quantization has 
resulted into a variety of scenarios in which a classical singularity is replaced by 
a quantum mechanical bounce, both in the cosmological context El [HI [91 [10] 
and with black holes [HJ H21 OSl M H5J US] • 

In this paper we investigate polymer quantization of spherically symmetric 
Einstein gravity in a description in which the classical dynamical variables are 
adapted to the evolution of the Einstein- Rosen wormhole throat [T7J [THJ [19] • As 
expected from Birkhoff 's theorem and from Hamiltonian analyses of spherically 
symmetric Einstein gravity [20| |2T| 122] , the phase space of the reduced Hamilto- 
nian theory is two-dimensional, and the configuration variable specifies the size of 
the wormhole throat as it evolves from the past singularity through the bifurcation 
two-sphere to the final singularity. (For generalizations to charged and rotating 
black holes, see [23].) The advantage of this model, compared with a number of 
recent black hole mini-superspace models based on a spatially homogeneous foli- 
ation [TTJ [12j [131 CHI CHS EB] ) is that the reduced Hamiltonian is sufficiently simple 
to allow not just a study of the effective classical polymer dynamics but also a 
numerical extraction of the mass eigenvalues. Further, the model does not rely on 
auxiliary structures that are typically encountered in mini-superspace models of 
cosmology and black holes [9]. A disadvantage of the model however is that the 
reconstruction of a full quantum-corrected spacetime would need further input, 
including a specification of where the spatial hypersurfaces of the Hamiltonian 
foliation intersect the horizons. 

One piece of input that the polymerized throat dynamics model does require 
is the choice of the phase space variable that will be polymerized. It has been 
observed within cosmological mini-superspace models that this choice can signif- 
icantly affect the polymer dynamics, even in its qualitative properties [61 [9] . We 
show that the choice of the polymerized variable has significant consequences for 
the throat dynamics already at the semiclassical level, and in particular only cer- 
tain choices replace the classical singularity by a quantum bounce in the effective 
polymer dynamics. We then focus on the case where the polymerized variable is 
the area of the wormhole throat. This choice is motivated by the proposal of an 
evenly-spaced area spectrum for quantum black holes, obtained from a wide range 
of premises; see [2U E3 ESI EZ1 ESI ESI ESI ED] for a selection, and the references 
in [19] for a more extensive list. We show that polymerizing the wormhole throat 
area does lead to a bounce. A perhaps unexpected property however is that for 
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a macroscopic black hole the scale of the bounce can be significantly above the 
Planck scale, unless the polymerization scale itself is taken to be well below the 
Planck scale. 

From a mathematical viewpoint, a qualitative difference between our system 
and the 1/r and 1/r 2 potentials considered in [21 H] is that our polymer Hamil- 
tonian is nonsingular, and there is no need to introduce a "Thiemann trick" [31] 
modification in the Hamiltonian to regularize any singular terms. The eigenstate 
recursion relation obtained from our Hamiltonian can contain a singularity, with 
certain factor orderings, but within the orderings that we consider, this singular- 
ity acts merely as a boundary condition in the recursion relation. Reassuringly, 
we find that the semiclassical regime of the theory is not sensitive to the factor 
ordering. 

Our main results come from a numerical evaluation of the spectrum. We show 
that the distribution of the large eigenvalues implies an evenly-spaced area spec- 
trum, and these large eigenvalues are in good agreement with a Bohr-Sommerfeld 
semiclassical estimate. The choice of the factor ordering has a significant effect 
only on the lowest few eigenvalues. In the limit of small polymerization scale 
we recover, within the numerical accuracy, the area spectrum obtained from a 
Schrodinger quantization of the wormhole throat dynamics. We have not at- 
tempted to support the numerical results by a rigorous investigation of the spec- 
tral properties of the polymer Hamiltonian, but comparison with the correspond- 
ing Schrodinger Hamiltonian analysis [19] strongly suggests that our polymer 
Hamiltonian is essentially self-adjoint and has a discrete spectrum. 

The paper is organized as follows. In Sec. [2] we recall the main features of 
the classical throat theory and introduce a canonical chart that is suitable for 
polymerization. A conventional Schrodinger quantization is performed in Sec. El 
In Sec. H] we give a brief review of polymer quantization, including the classical 
limit known as effective polymerization, and the effective polymerization of the 
throat theory is discussed in Sec. |5j Full polymer quantum theory is analysed 
in Sec. El including the numerical results for the mass eigenvalues. Section [7] 
presents a summary and brief concluding remarks. 

Unless otherwise stated, we use Planck units in which G — c — h — 1. 



We start with a Hamiltonian system with a two-dimensional phase space and the 
Hamiltonian 



where the configuration variable r takes positive values and p is the conjugate 
momentum. The equations of motion reduce to 



2 Classical throat theory 





2M 



(2.2) 
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where M is the conserved value of H and the overdot denotes derivative with 
respect to the time t. As observed in [T71 [18], (I2.2p is the equation of a radial 
timelike geodesic that passes through the bifurcation two-sphere on a Kruskal 
manifold of mass M, with r being the area-radius of the two-sphere and t the 
proper time. It was shown in [TH] that the Hamiltonian ( 12.1 p can be obtained by 
a Hamiltonian reduction of the spherically symmetric sector of Einstein's theory 
under suitable boundary conditions: the spacelike hypersurfaces are frozen at 
the (say) left-hand-side i° of the Kruskal diagram, they evolve at unit rate with 
respect to the asymptotic Minkowski time at the right- hand- side i°, and they 
intersect a radial geodesic through the bifurcation two-sphere so that t coincides 
with the proper time on this geodesic. The variable r can then be identified as 
area-radius of the two-sphere on the distinguished geodesic. If in addition the 
spatial slices are chosen to intersect the radial timelike geodesic at the point where 
the two-sphere radius is minimized, the Hamiltonian (12. ip can be regarded as the 
Hamiltonian of the Einstein- Rosen wormhole throat. 

The Hamiltonian system thus describes the proper time evolution of the 
Einstein-Rosen wormhole throat radius as it expands from zero at the past sin- 
gularity, reaches the maximum value 2M at the bifurcation two- sphere, and col- 
lapses back to zero at the future singularity. This means that the spacetime 
dynamics is described in terms of variables that are in a certain sense confined 
'inside' the black hole. This is qualitatively similar to the recent studies of 
quantum black holes in terms of spatially homogeneous slicings of the interior 
[TT| [T2| IT3j [TJJ [151 EE], but there are two significant differences. First, our sys- 
tem is fully reduced, with a true Hamiltonian and no constraints. This is an 
advantage in the sense that implementation of polymer quantization will be rel- 
atively straightforward, and it will be easy to compare the results to those of 
a conventional Schrodinger quantization. Second, the foliation of our system is 
specified at the infinities and at the wormhole throat, but the foliation is largely 
arbitrary in the intermediate regions, including the locations where the spacelike 
hypersurfaces cross the future and past branches of the Killing horizon. This arbi- 
trariness is a disadvantage in the sense that additional input on the foliation seems 
necessary before one could attempt to reconstruct an entire quantum-corrected 
spacetime from our quantum theory. 

The purpose of this paper is to examine polymer quantization of the wormhole 
throat Hamiltonian ( 12. ip . As the results of polymer quantization can depend 
significantly on the precise choice of the polymerized variables [S], we shall 
set up the problem in a way that allows polymerization of an arbitrary positive 
power of the radius: we introduce the new canonical chart ((f), II), where 




(2.3a) 
(2.3b) 



n = apr 



l-l/a 
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and a is a positive constant. The Hamiltonian takes the form 



h = t{^ + 1 )- < 2 - 4 > 

The case a = 1/2 is geometrically special in that </> is then proportional to the 
throat area. This case is also special in that if we had retained geometric units, 
in which G = c = 1 but [H] = (length) 2 , II would be dimensionless precisely for 
a = 1/2. 



3 Schrodinger quantization 

As the classical phase space is finite dimensional, the Hilbert space of a conven- 
tional Schrodinger quantization can be chosen to be square integrable functions 
on a finite-dimensional configuration space. We choose the configuration variable 
to be 0, and we take the inner product to have a flat measure in 0: 



W>iM= / M<l>)M<!>)d<j>. (3.1) 



o 

To quantize the Hamiltonian (12. 4p . we make the replacement II i— > —id/dcf) and 
choose a symmetric factor ordering. For reasons that will emerge in section [6l we 
consider the family of symmetric orderings in which 

where the ordering parameter (3 takes real values. 

Writing = x 2 /( 3a ) and = x^ 1 ^ 2 " > ~ l ^' io ^x{ x )-, "this Schrodinger theory is 
mapped to one in which the inner product is 



(Xi,X2)o= Xi(x)X2(x)dx, (3.3) 



o 



where we have omitted an overall multiplicative constant, and the Hamiltonian 
is 

9 / d 2 Ax 2 /* (9a + 4/3-2) (3a + 4/3-2) 



ffo l["tf + T + 36W ) ■ (3 ' 4) 

This is the theory discussed in [19], and the parameter denoted in [19] by r has 
the value (9a + 4/3 — 2) /(6a). H is essentially self-adjoint for (3 > ■= and for 
(3 < \ — 3a, and for other values of of the parameters Ho has a U(l) family of 
self-adjoint extensions, characterized by a boundary condition at x — 0. In all 
cases the spectrum of Hq is discrete, and the large eigenenergies have the WKB 
estimate 

E% fKB ~2k + B + o(\), (3.5) 
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where k is an integer and o(l) denotes a term that vanishes asymptotically at 
large E. The additive constant B depends on a and (3 and, for -= — 3a < /3 < tt, 
also on the self-adjoint extension. For certain choices of the self-adjoint extension, 
the WKB estimate differs from the true eigenvalues by only a few percent already 
near the ground state [32J. 



4 Polymer quantization: an overview 

Polymer quantization is a loop quantum gravity motivated, unitarily inequiva- 
lent alternative to conventional Schrodinger quantization [TJ [2] . It incorporates 
fundamental discreteness to the underlying spatial geometry by confining the 
quantized Hamiltonian dynamics into a discrete spatial lattice. While the overall 
philosophy of polymer quantization is deeply connected to loop quantum gravity, 
polymer quantization is a consistent quantization scheme in its own right, and it 
is of interest apply this scheme to a variety of systems. 

Consider a quantum theory on the real line R. The polymer Hilbert space 
is spanned by the normalizable eigenstates |x) of the position operator with the 
inner product 

(x'|x) = 5 x , x , (4.1) 

where the quantity on the right hand side is the Kronecker delta. Any state in 
the polymer Hilbert space has support on only countably many points, which 
makes it impossible to define the momentum operator as a differential operator. 
Instead one defines the actions of two basic operators, the position operator x 
and the finite translation operator U^, as 

x |x) = x |x), (4.2a) 
U,\x) = \x + fj), (4.2b) 

where \x > is the polymerization scale. In general \i can be a function of x but 
in all that follows it is assumed to be a constant. The momentum operator is 
then constructed from the translation operator [I]: 

#=^-e»- (4.3) 

In the \i — > limit, (14. 3D reduces to the standard momentum operator p = 
—id x , thus leading to the usual Schrodinger quantized system [33]. For critical 
discussions of polymer quantum mechanics, see [34] . 

The (semi) classical regime of polymer quantum theory is obtained by keeping 
the polymerization scale \i fixed and making the replacement t— > e w , where 
p is the continuum momentum. This means making in the classical continuum 
Hamiltonian the replacement p h- > sin(yup) / p,. This limit, known as the effective 
polymer theory or effective polymerization, has proved a useful tool with mini- 
superspace models of quantum cosmology and black holes 0, 0, El El [HI [12l [T3"| 

mug aging. 
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5 Effective polymer theory 



We now turn to the effective polymerization of the Hamiltonian ( 12.41) . We shall 
see that the effective theory will restrict the physically interesting values of the 
parameter a. 

The effective Hamiltonian is obtained from (12. 4p by the replacement II y 
sin(yull)/yu: 

^ / ^ 2 - 4a sin 2 ( yU n) 
2 v «V 

H takes positive values and is a constant of motion. Denoting the value of H 
by M, we see that is bounded above by < (2M) 1//a =: + . Note that 
this upper bound is at the turning point of the non-polymerized motion, at the 
bifurcate two-sphere of the black hole. 

When a > 2/3, there are solutions that reach = 0. We shall not consider 
this case further. 

When a < 2/3, the motion is oscillatory, with the outer turning point at 
= 0+ and the inner turning point at = 0_, where 0_ is the unique solution 
to 

A2-3a 

2M = ^-- + a . (5.2) 

The outer turning point is at the turning point of the non-polymerized motion, but 
the inner turning point is a genuine polymerization effect, replacing the singularity 
of the non-polymerized theory by a bounce. There are three qualitatively different 
subcases: 

1. l/2<a<2/3. The asymptotic forms of 0_ at large and small M are 

0_ « (2M) 1/a , M -> oo, (5.3a) 

0_ w (2M« 2 /i 2 ) 1/(2 ~ 3a) , M -» 0. (5.3b) 



2. a = 1/2. The inner turning point is 

2 \ 2 



,4 + // 2 

3. a < 1/2. The asymptotic forms of 0_ at large and small M are 



(5.4) 



0_ w (2MaV) 1/(2 ~ 3a) , M — ^ oo, (5.5a) 

0_ « (2M) 1/a , M^O. (5.5b) 

In case 1 we have 0_/0 + — >• 1 as M — > oo. This means that for a macroscopic 
black hole the bounce happens close to the outer turning point, which seems 
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physically undesirable. In case 2 the ratio 0_/0 + is independent of M, while in 
case 3 we have 0_ — > oo but 0_/0+ — > as M — > oo. Cases 2 and 3 seem hence 
more reasonable, and taking fi a few orders below unity will in these cases result 
into a bounce radius much smaller than the size of the black hole. It is perhaps 
disconcerting that if /i is of order unity, the bounce in cases 2 and 3 may occur 
at a macroscopic scale, and in case 2 even at a scale that is comparable to that 
of the hole. However, macroscopic polymerization effects have been encountered 
also in the cosmological context El El El [10] , and in the absence of a way 
to recover from the throat theory a full quantum-corrected spacetime, it is not 
obvious that these properties would be physically unacceptable. We therefore 
regard cases 2 and 3 as physically interesting. 

In the rest of the paper we will focus on a = 1/2 for several reasons. First of 
all, polymerizing the area is motivated by the number of quantization methods 
that lead to an evenly-spaced area spectrum of a quantum black hole [2H [251 [261 
|2"7] dHl ESI EHl ED], including Schrodinger quantization of the throat theory [15] , 
as seen in f)3.5p . One therefore expects the qualitative similarities and differences 
between polymer and Schrodinger spectra to be most easily detectable for a = 
1/2. We speculate that for a < 1/2 the polymer mass spectrum may differ more 
drastically from the Schrodinger spectrum, although we still expect the two to 
agree in the limit fi — > 0. 

Second of all, with a = 1/2 we can compare the numerical results to an 
analytic evaluation of the large eigenvalues by the Bohr- Sommerf eld semiclassical 
quantization rule. The Bohr- Sommerf eld rule states that adiabatic invariants of a 
classical system in periodic motion are quantized in units of Planck's constant h. 
In our system the canonical variable with periodic time evolution is <fr, whereas IT 
is monotonically decreasing, and the turning points occur when II is an integer 
multiple of 7r/(2/i). The relevant adiabatic invariant is therefore — § 011 dt, where 
the integral is evaluated over one period in 0. Setting a = 1/2, we have 



where n is a non-negative integer, we have written /ill = 6, and in the last equality 
we have evaluated the elementary integral. Recalling that in our units h = 2tt, 
(15.61) gives for M 2 the spectrum 



nh 





(5.6) 




3/2 



M 2 = k{jj)n, 



(5.7a) 



fc(Ai): 




(5.7b) 



2(2 + /i 2 ) ' 
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The key observation about (j5.7p is it that it implies an evenly-spaced area 
spectrum, with the /i-dependent spacing factor k(fi). In the continuum limit 
jj, — y we have k(fi) — > 2, and the spectrum agrees with the leading term in the 
Schrodinger spectrum (13.51) . The spacing factor has a minimum at /i = a/2, and 
for large \i it asymptotes to /i/2. We will see in section IfTBI that this spectrum is in 
excellent agreement with our numerical results for the full polymer quantization 
except for the lowest few eigenvalues. 



6 Full quantum polymer theory 

We now proceed to the fully polymerized quantum throat theory. The main aim 
will be a numerical evaluation of the spectrum. 

For the reasons discussed in Sec. |5j we set a = 1/2, so that the polymerized 
variable is proportional to the throat area. The quantum polymer Hamiltonian 
is obtained from (12. 4p as described in Sec. HJ acts on the eigenstates of by 
multiplication, while 

n:=^(^-^), (6.1a) 

U lt \(f>) = \<f> + fi), (6.1b) 

where fi > is the same polymerization scale that appeared in the effective 
polymer theory of Sec. |5j In order to deal with negative values of and fractional 
powers of 0, we define 

fr\4>) := |0p|0> (6.2) 

for any real number 7. 

The kinetic term in the Hamiltonian presents an issue of factor ordering. We 
consider two orderings in turn. 



6.1 Nonsingular factor ordering 

In this subsection we order the polymer Hamiltonian as 

H pol = ~ (ffl^n + 1/2 ), (6-3) 

corresponding to the Schrodinger Hamiltonian (13. 2 p with a = 1/2 and /3 = 0. For 
reasons that will emerge shortly, we refer to (16. 3p as the nonsingular ordering. 

We write the basis states as |0) := \mjl), where m e Z and Jx := 2/i. We then 
have 

-2 



if po i I m/i) 



3/2 



m + I 



1/2 



\m 



1 1 1/2 
2 1 



4 



m 



1/2 



\ m + 2 



11 1/2 



\{m + l)p) 



m- 2 



11 1/2 



|m/i) 



|(m - l)/i) 



(6.4) 



8 



The subspace in which c m = c_ m and the subspace in which c m = — c_ m are 
each invariant under the action of the Hamiltonian. These subspaces therefore 
decouple and it suffices to consider each of them individually. We refer to these 
subspaces respectively as the symmetric sector and the antisymmetric sector. 
The eigenvalue equation 



H po \ 22, c ™\ m v) = M c m \mfi) 

m=— oo m=— oo 

rives the recursion relation 



(6.5) 



m + 



II 1/2 



+ m- 



1 1 1/2 

2 I 



4 



\m 



1/2 _ 



\m- 



i\V 2 , 



_! + |m+2| c m+1 . 

(6.6) 

As the coefficients of c m _i and c m+ \ in (16. 6 p are nonvanishing for all m, the 
recursion relation is nonsingular in the sense that the free initial data consist 
of c m at any two adjacent values of m, and these initial data determine c m for 
all m. Expecting a discrete spectrum, we look for solutions that are normalizable, 



Er 



< oo. Examination of the asymptotic form of (16. 6p shows that such 



solutions have at m — > ±oo the asymptotic form 



c m = \m\ 5 1 ^ exp(A^/|?ri|) 



i + ^ + ,/ri + ^ 2 -i 



l + 0(|mr 1/2 ) 



(6.7) 



where 



d = x 

32 



1 + 



2\ 2 



3/2 ' 



A 



3/2 M 



i + r 2 



1/2" 



(6.8a) 



(6.8b) 



To find the mass eigenvalues numerically we use the shooting method. Given 
a trial value for M, we choose an mo so large that c mo and c mo _i can be obtained 
from (16 .71) and iterate downwards using (16. 61) . In the symmetric sector we shoot 
for values of M that give c\ = c_i, and in the antisymmetric sector we shoot for 
Co = 0. The numerical accuracy is monitored by increasing mo until the results 
no longer change to the desired accuracy. 

Working to six significant figures, we find reliably the 10 lowest eigenvalues 
for jj, < 15 and the 50 lowest eigenvalues for fi < 2.5. For small /i the numerics 
becomes increasingly slow and we do not consider values of \i below 0.005. 
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0.5 



1.5 



2.5 



M 



Figure 1: Five lowest symmetric sector mass eigenvalues for 0.005 < fi < 2.5. 
(The sizes of the markers do not indicate the numerical uncertainty.) 




Figure 2: Five lowest symmetric sector mass eigenvalues for 0.025 < \x < 10. 
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Figures [T] and [2] show the lowest five eigenvalues for 0.005 < \i < 10 in the 
symmetric sector, and figures [3] and H] show the corresponding data in the anti- 
symmetric sector. The ground state is in the symmetric sector, as expected, and 
the symmetric and antisymmetric sector eigenvalues alternate. The ground state 
eigenvalue decreases with increasing fi and appears to asymptote to zero from 
above as fi — > oo. The higher eigenvalues show a more complicated dependence 
on /j, for /j, < 2, most of them having a minimum that is close to the minimum of 
the function k(fi) (15.7b|) . but they all increase in \x for \x > 2, at a rate that we 



will examine in Sec. 16.31 The symmetric and antisymmetric sector eigenvalues 
do not appear to asymptote to each other as fi — > oo, but for fixed fi the gap 
between the symmetric and antisymmetric sector eigenvalues gets smaller as the 
mass level increases. 

In the \i — > limit the eigenvalues appear to converge to definite limiting 
values. We wish to compare these limiting values to the eigenvalues in the corre- 
sponding Schrodinger theory, with the inner product (13.11) and the Hamiltonian 
(13. 2p with a = 1/2 and /3 = 0. Mapping this Schrodinger theory to that of (13. 3p 
and (13 .4p , discussed in [19] , we find that the parameter denoted in [19] by r has the 
value 5/6. The eigenf unctions are parabolic cylinder functions and the eigenvalue 
equation is (B10) of [19]. The small argument behaviour of the eigenf unctions 
is = cos(^)0 1 / 2 [1 + O(0 3/2 )] + sin(0) [1 + O(0 3 / 2 )], where 8 G [0,tt) is the 
parameter that specifies the self-adjoint extension of the Hamiltonian. This small 
argument behaviour suggests that the symmetric and antisymmetric sector of the 
polymer theory should correspond respectively to the 8 = n/2 and 8 = Hamil- 
tonians in the Schrodinger theory. With this identification of 8, we have verified 
that the polymer theory eigenvalues indeed converge to those of the correspond- 
ing Schrodinger theory, for the low eigenvalues where numerical accuracy allows 
us to test the convergence. The evidence for the three lowest eigenvalues in each 
sector is shown in Table [U 

6.2 Factor ordering with a singular eigenvalue equation 

In this subsection we order the polymer Hamiltonian as 

H pol =l(4V*flH^ + ^), (6.9) 

corresponding to the Schrodinger Hamiltonian (13. 2p with a = 1/2 and ft = 1/4. 
For reasons that are about to emerge, we refer to ( 16. 9p as the singular ordering. 
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Nonsingular ordering, symmetric sector 





M 


M l 


M 2 


0.1 


0.763255 


1.62108 


2.15328 


0.05 


0.766485 


1.62254 


2.15461 


0.005 


0.768132 


1.62321 


2.15517 


9 = 7r/2 Schrodinger 


0.768184 


1.62323 


2.15518 



Nonsingular ordering, antisymmetric sector 



n 


M 




Mi 


M 2 


0.1 


1.19141 


1. 


85036 


2.32984 


0.05 


1.18178 


1. 


S4343 


2.32411 


0.005 


1.16602 


1. 


83184 


2.31424 


9 = Schrodinger 


1.15890 


1. 


82661 


2.30978 



Singular ordering 



n 


M 


Mi 


M 2 


0.1 


1.31307 


1.92381 


2.38555 


0.05 


1.31320 


1.92415 


2.38605 


0.005 


1.31325 


1.92427 


2.38623 



Table 1: Three lowest eigenvalues in the symmetric and antisymmetric sec- 
tors with the nonsingular ordering and with the singular ordering as \i — > 0. 
For the nonsingular ordering, the eigenvalues convergence to the corresponding 
Schrodinger eigenvalues. 



Writing again /x = 2/x, we now have 
2 



if po i|m/i) 



3/2 



2 + — ) \m\ \mfi) 



\m{m + l)| 1//4 |(m + l)/x) — \m(m — lJj^Km — l)/i) 



(6.10) 

The action of the Hamiltonian leaves the subspaces of positive m and negative 
m invariant and annihilates the subspace m = 0. These three subspaces are thus 
entirely decoupled. The dynamics in the m = subspace is trivial, and the only 
energy eigenvalue is zero. The actions in the subspaces of positive and negative 
m are isomorphic via m t— > —m, and it hence suffices to consider the subspace of 
positive m. 
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The eigenvalue equation 



oo oo 



(6.11) 



m=l m=l 



gives the recursion relation 




2 



= m 



■m—l 



(6.12) 



where m — 1,2,3, ... . Note that this relation does not involve cq as its coefficient 
vanishes. The recursion relation is singular in the sense that evaluating ( 16 . 1 2 j) for 
m = 1 gives the condition 



and the free initial data hence consists of only c±. Examination of the asymptotic 
form of (I6.12p at m — > oo shows that the normalizable solutions have at m — > oo 
the asymptotic behaviour given by (16. 7p and (16. 8p . 

The numerical problem is almost identical to that of Sec. 16.11 We use (I6.12p to 
come down from a large m , matching c mo and c mo _i to (16. 7p . and shoot for ( I6.13p . 

Figures [5] and [6] show the lowest five eigenvalues for 0.005 < \i < 10. Compar- 
ison with Figs. |3] and H] shows that the eigenvalues are qualitatively very similar 
to those in the antisymmetric sector of the nonsingular ordering, albeit slightly 
higher: the two spectra appear to asymptote to each other as ji — > oo, and the 
convergence gets more rapid as the mass level increases. At /i — > 0, the eigen- 
values appear to have well-defined limits that are genuinely above those of the 
antisymmetric sector of the nonsingular ordering, as illustrated in Fig. [7] and in 
Table [TJ We expect this limit to give the eigenvalues of the Schrodinger Hamil- 
tonian (13. 2 p with a = 1/2 and (3 = 1/4, with an appropriate self-adjointness 
boundary condition, although we have not attempted to verify this. 

In summary, the spectrum of the singular ordering Hamiltonian (I6.9P is quali- 
tatively very similar to the spectrum of the antisymmetric sector of the nonsingu- 
lar ordering Hamiltonian (16. 3p . The boundary condition (I6.13P that is enforced by 
the singular recursion relation (I6.12p plays the part of restricting the nonsingular 
recursion relation (16. 6 p to the antisymmetric sector. 

6.3 Semiclassical spectrum 

We now turn to the semiclassical regime of the spectra. The Schrodinger quanti- 
zation asymptotic estimate (I3.5P and the Bohr- Sommerf eld polymer result ( 15. 7p 
suggest that the asymptotic behaviour of the large eigenvalues is 




(6.13) 



Ml ~ k(ji)n + B, 



(6.14) 
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Figure 7: Two lowest eigenvalues for the antisymmetric sector of the nonsingular 
ordering (purple curve with square markers) and for the singular ordering (dark 
blue curve with diamond markers). 



5 




1 A 1 r- 1 1 

2 4 6 8 

M 

Figure 8: The spacing factor k(n) as obtained from the numerical eigenvalues (dis- 
crete points) and from the Bohr- Sommerf eld estimate (15.71) (continuous curve). 
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where k(fi) is the Bohr-Sommerfeld spacing function (15.7bl) . for both sectors of 
the nonsingular ordering and for the singular ordering. The additive constant B 
could a priori depend on the ordering, the sector and \i. 

Our numerical results confirm this expectation in the range of fi where the 
numerical accuracy is sufficiently good. The numerical limitations arise mainly 
from the computational noise that builds up after 50 or so eigenvalues, or for 
fi ~ 15 already after 15 or so eigenvalues. For 0.65 < fi < 15 we find full 
agreement with (I6.14p . For fi < 0.65 the noise piles up before the semiclassical 
regime is reached, with the notable exception of fi ~ 0.01, where the spectrum 
appears to obey f)6.14p already near the ground state with k very close to 2. 
A plot of the numerical results superposed on the Bohr-Sommerfeld k(fi) is shown 
in Fig. E 

The numerical results give also information on the additive constant B for 
0.65 < \i < 15. With the nonsingular ordering, we find that B differs between 
the two sectors. In the symmetric sector B is small and negative (taking the 
convention in which the ground state is n = 0), and it is essentially independent 
of \i. In the antisymmetric sector B is positive, and it depends on /i qualitatively 
in the same way as k(fi). With the singular ordering, B is very close to that 
of the nonsingular ordering antisymmetric sector, and the difference appears to 
be decreasing as fi increases (cf. the behaviour of the lowest two eigenvalues in 
Fig. [7]). The difference is however small, and we have not been able to exclude 
the possibility that it might be an artefact of the numerics. 

7 Conclusions 

In this paper we have addressed polymer quantization of spherically symmetric 
Einstein gravity, in a description in which the dynamical variables are associ- 
ated with the Einstein- Rosen wormhole throat. Choosing to polymerize the area 
of the throat, we found numerically that the spectrum is discrete and positive 
definite, and the large eigenvalues imply an evenly-spaced area spectrum that is 
consistent with a Bohr-Sommerfeld semiclassical estimate. Issues of factor or- 
dering arise in the kinetic term of the Hamiltonian, but the spectrum is largely 
insensitive to these issues except at the lowest few eigenvalues. In the limit of 
small polymerization scale the spectrum tends, within the numerical accuracy, to 
that of a conventional Schrodinger quantization. 

Our mathematical results are broadly similar to those obtained previously in 
polymer quantizations of the 1/r and 1/r 2 potentials [31 H], despite the technical 
differences that arise from the singularities of these potentials. There is however a 
difference in the physical interpretation of the polymerization scale. The physical 
context of the 1/r potential is quantum gravity corrections to atomic physics, and 
in this case we expect the polymerization scale fi to be at the Planck scale, which 
is much smaller than the Rydberg scale of the non-polymerized Hamiltonian. It 
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was indeed found in [3] that the polymer results are in this limit close to the 
usual Schrodinger results. The physical context of the 1/r 2 potential is less clear, 
but as this system is classically scale invariant, the polymerization scale \x affects 
the results only by setting the energy scale. Now, the physical context of our 
system is genuine gravitational dynamics, and one expects the polymerization 
scale fi to be at the Planck scale, which is of order unity in our Planck units. Our 
numerical results for the spectrum appear reasonable for /i ~ 1, but it is puzzling 
that the wormhole throat bounce in the effective polymer theory then occurs at 
a macroscopic scale, for a macroscopic black hole. It would be of interest to 
examine the motion of wave packets in the full polymer theory to ascertain that 
this macroscopic bounce is not an artefact of the effective polymer theory. 

The bounce necessarily occurs below the horizon, so the question of whether a 
macroscopic bounce scale has observational consequences for a galactic black hole, 
for example, depends on the details of the quantum corrected space-time exterior 
to the horizon. Reconstructing a quantum-corrected spacetime from the polymer 
throat theory would require additional input about the spatial hypersurfaces of 
the Hamiltonian foliation. It might be of interest to analyze this reconstruction 
in some geometrically simple foliations, such as the Novikov coordinates [35], 
but the absence of local dynamical degrees of freedom would probably make it 
difficult to prefer any such foliation to another. The prospects of fixing a foliation 
by local criteria would be higher if our techniques could be generalized to systems 
that carry Hawking radiation, such as the spherically symmetric Einstein-scalar 
system [36J. In the absence of such additional structure, however, one would 
hope that local changes in the foliation (at least those that preserve spherical 
symmetry) would not affect the physical predictions of the model. If this were 
not the case, it would reinstate in a much more serious way the problem of 
dependence on fiducial structures that occurs in the semi-classical polymerization 
of the homogeneous interior. 

Acknowledgments 

We thank Jonathan Ziprick for helpful discussions. GK was supported in part 
by the Natural Sciences and Engineering Research Council of Canada. JL was 
supported in part by STFC (UK) grant PP/D507358/1. 

References 

[1] A. Ashtekar, S. Fairhurst and J. Willis, "Quantum Gravity, Shadow 
States, and Quantum Mechanics", Class. Quant. Grav. 20, 1031 (2003) 
|arXiv:gr-qc/0207106] . 



18 



[2] H. Halvorson, "Complementarity of Representations in Quan- 
tum Mechanics", Studi es Hist. Philos. Mod. Phys. 35, 45 (2004) 
|arXiv:quant-ph/0110l02] . 

[3] V. Husain, J. Louko and O. Winkler, "Quantum Gravity and the Coulomb 
Potential", Phys. Rev. D 76, 084002 (2007) |arXiv:0707.0273l [gr-qc]]. 

[4] G. Kunstatter, J. Louko and J. Ziprick, "Polymer Quantization, Singular- 
ity Resolution and the 1/r 2 Potential", Phys. Rev. A 79, 032104 (2009) 
|arXiv:0809.5"098l [gr-qc]]. 

[5] D. Green and W. G. Unruh, "Difficulties with Recollapsing models in Closed 
Isotropic Loop Quan tum Cosmology", Phys. Rev. D 70, 103502 (2004) 
|arXiv:gr-qc/0408074"] . 

[6] A. Ashtekar, T. Pawlowski and P. Singh, "Quantum Nature of the 
Big Bang: Improve d Dynamics", Phys. Rev. D 74, 084003 (2006) 
|arXiv:gr-qc/0607039l ; 

[7] A. Ashtekar, T. Pawlowski, P. Singh and K. Vandersloot, "Loop Quan- 
tum Cosmology of k = 1 FRW Models", Phys. Rev. D 75, 024035 (2007) 
arXiv:gr-qc/0612104"] ; 



[8] K. Vandersloot, "Loop Quantum Cosmology and the k = -1 RW Model", 
Phys. Rev. D 75, 023523 (2007) |arXiv:gr-qc/0612070] . 

[9] A. Corichi and P. Singh, "Is Loop Quantization in Cosmology Unique?", 
Phys. Rev. D 78, 024034 (2008) |arXiv:0805.0~l36 j: "A Geometric Perspective 
on Singularity Resolution and Uniqueness in Loop Quantum Cosmology", 
Phys. Rev. D 80, 044024 (2009) jarXiv: 0905 .49491 [gr-qc]]. 

[10] Y. Ding, Y. Ma and J. Yang, "Effective Scenario of Loop Quantum Cosmol- 
ogy", Phys. Rev. Lett. 102, 051301 (2009) [ arXiv:0808.0990l [gr-qc]]. 

[11] A. Ashtekar and M. Bojowald, "Quantum Geometry and the Schwarzschild 
Singularity", Class. Quant. Grav. 23, 391 (2006), |arXiv:gr-qc/0509075] . 

[12] L. Modesto, "Loop Quantum Black Hole", Class. Quant. Grav. 23, 
5587 (2006) [arXiv:gr-qc/0509078] "Black Hole Interior from Loop 
Quantum Gravity", Adv. High Energy Phys. 2008, 459290 (2008) 
[arXiv:gr-qc/0611043| ; "Space- Time Structure of Loop Quantum Black 
Hole" , larXiv:0811.2T96l [gr-qc]. 

[13] C. G. Boehmer and K. Vandersloot, "Loop Quantum Dynamics of the 
Schwarzschild Interior", Phys. Rev. D 76, 104030 (2007) [ arXiv:0709.2T29l 
[gr-qc]]; "Stability of the Schwarzschild Interior in Loop Quantum Gravity", 
Phys. Rev. D 78, 067501 f2008) larXiv":0807.3042l [gr-qc]. 



19 



[14] M. Campiglia, R. Gambini and J. Pullin, "Loop Quantization of Spherically 
Symmetric Midi-Superspaces : The Interior Problem" , AIP Conf. Proc. 977, 
52 (2008) [arXiv:0712.0817l [gr-qc]]: R. Gambini and J. Pullin, "Black holes 
in Loop Quantum Gravity: the Complete Space-Time", Phys. Rev. Lett. 
101, 161301 (2008) [arXiv: 0805. 1187] . 

[15] W. Nelson and M. Sakellariadou, "Numerical Techniques for Solving the 
Quantum Constraint Equation of Generic Lattice-Refined Models in Loop 
Quantum Cosmology", Phys. Rev. D 78, 024030 (2008) |arXiv:0803.4483l 
[gr-qc]]. 

[16] A. Peltola and G. Kunstatter, "A Complete, Single-Horizon Quantum 
Corrected Black Hole Spacetime", Phys. Rev. D 79, 061501(R) (2009) 
|arXiv: 08 1 1 . 3240 [gr-qc]]; "Effective Polymer Dynamics of D-Dimensional 
Black Hole Interiors", Phys. Rev. D 80, 044031 (2009) |arXiv:0902. 17461 
[gr-qc]]. 

[17] J. L. Friedman, I. H. Redmount, and S. N. Winters-Hilt (unpublished). 

[18] I. H. Redmount, talk given at the 3rd Midwest Relativity Meeting, Oakland, 
Michigan, 1993. 

[19] J. Louko and J. Makela, "Area spectrum of the Schwarzschild black hole", 
Phys. Rev. D 54, 4982 (1996) |arXiv:gr-qc/9605058] . 

[20] T. Thiemann and H. A. Kastrup, "Canonical quantization of spherically 
symmetric gravity in Ashtekar's selfdual representation," Nucl. Phys. B 399, 
211 (1993) |arXiv:gr-qc/9310012] . 

[21] H. A. Kastrup and T. Thiemann, "Spherically Symmetric Gravity As 
A Completely Integrable System," Nucl. Phys. B 425, 665 (1994) 
|arXiv:gr-qc/9401032] . 

[22] K. V. Kuchar, "Geometrodynamics of Schwarzschild black holes," Phys. Rev. 
D 50, 3961 (1994) |arXiv:gr-qc/9403003] . 

[23] J. Makela and P. Repo, "A Quantum Mechanical Model of the Reissner- 
Nordstrdm Black Hole", Phys. Rev. D 57, 4899 (1998) |arXiv:gr-qc/9708029] ; 
J. Makela, P. Repo, M. Luomajoki and J. Piilonen, "Quantum-Mechanical 
Model of the Kerr-Ne wman Black Hole", Phys. Rev. D 64, 024018 (2001) 
|arXiv:gr-qc/0012055] . 

[24] J. Bekenstein, "The Quantum Mass Spectrum of the Kerr Black Hole", Lett. 
Nuovo Cim. 11, 467 (1974). 

[25] J. D. Bekenstein and V. F. Mukhanov, "Spectroscopy of the Quantum Black 
Hole", Phys. Lett. B 360,7 (1995) |arXiv:gr-qc/9505012] . 



20 



H. A. Kastrup, "On the Quantum Levels of Isolated Spherically Symmetric 
Gravitational Systems", Phys. Lett. B 385, 75 (1996) |arXiv:gr-qc/9605038] . 

A. Barvinsky and G. Kunstatter, "Exact Physical Black Hole States 
i n Generic 2-D Dila ton Gravity", Phys. Lett. B 389, 231 (1996) 
arXiv:hep-th/9606134] . 



A. Barvinsky, S. Das, G. Kunstatter, "Discrete Spectra of Charged Black 
Holes" , Found. Phys. 32, 1851 (2002) |arXiv:hep-th/0209039] 

G. Gour, A. J. M. Medved, "Quantum Spectrum for a Kerr-Newman Black 
Hole", Class. Quant. Grav. 20, 1661 (2003) |arXiv:gr-qc/0212021] . 

S. Hod, "Bohr's Correspondence Principle and The Area Spectrum of Quan- 
tum Black Holes", Phys. Rev. Lett. 81, 4293 (1998) |arXiv:gr-qc/9812002] . 

T. Thiemann, "Quantum spin dynamics (QSD)," Class. Quant. Grav. 15, 
839 (1998) [arXiv:gr-qc/9606089] . 

J. Makela and A. Peltola, "Spacetime Foam Model of the Schwarzschild 
Horizon", Phys. Rev. D 69, 124008 (2004) |arXiv:gr-qc/0307025] . 

A. Corichi, T. Vukasinac and J. A. Zapata, "Polymer Quantum Mechanics 
and Its Continuum Limit", Phys. Rev. D 76, 044016 (2007) jarXiv: 0704 .00071 
[gr-qc]]. 

H. Nicolai, K. Peeters and M. Zamaklar, "Loop Quantum Gravity: an Out- 
side View", Class. Quant. Grav. 22, R193 (2005) |arXiv:hep-th/0501114] ; 

R. C. Helling and G. Policastro, "String Quantization: Fock vs. LQG Rep- 



resentations" , arXiv:hep-th/0409182 



C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravitation (Freeman, San 
Francisco, 1973). 

J. Ziprick and G. Kunstatter, "Dynamical Singularity Resolution in Spher- 
ically Symmetric Black Hole Formation," Phys. Rev. D 80, 024032 (2009) 
|arXiv:0902.322l [gr-qc]]. 



21 



